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THE  POSSIBILITY  OP  USING  THE  TRANSPARENCY  METHOD 
TO  DETERMINE  THE  MICROSTRUCTURE 
OP  AN  ATMOSPHERIC  AEROSOL 

V.  G.  Bakhtiyarov 

Difficulties  in  Using  the  Transparency  Method 
in  the  Atmosphere 

The  transparency  method,  which  makes  it  possible  to  calculate  the 
spectra  of  the  particles  of  a  dispersed  system  according  to  optical 
information,  has  been  used  earlier  only  in  laboratory  conditions.  In 
this  work  we  consider  the  possibilities  of  using  the  method  in  the 
free  atmosphere. 

Recently  the  study  of  the  atmospheric  aerosol,  i.e.,  particles  of 
microscopic  and  submicroscopic  size  in  the  atmosphere,  has  taken  on 
ever  greater  importance.  In  connection  with  this  it  Is  necessary  to 
examine  the  question  of  the  possibility  of  applying  the  studied  method 
to  determine  the  dimensions  of  particles  of  an  atmospheric  aerosol. 

In  order  to  determine  particle  dimensions  by  the  transparency 
method  it  is  necessary  to  find  the  polydispersion  coefficient  of 
scattering  of  the  system  in  a  significant  range  of  wavelengths,  i.e., 
to  determine  the  transparency  of  the  atmosphere  in  a  certain  spectral 
interval . 

The  accepted  view  of  the  atmosphere  is  that  of  a  medium  consisting 
of  a  mixture  of  gases,  water  vapor,  and  extremely  fine  dispersed 
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particles  -  aerosols.  The  general  attenuation  of  radiation  in  the 
atmosphere  occurs  through  two  basic  processes:  absorption  by  gaseous 
components  and  scattering. 

Thus  the  attenuation  of  radiation  in  the  atmosphere  can  be  presented 
in  the  form 

QM-tM+«W.  (1) 

where  a(X)  is  the  attenuation  due  to  scattering  and  y(X)  is  the 
fundamental  absorption  of  the  atmosphere. 

The  dependence  of  y  on  wavelength  X  represents  a  series  of  lines, 
bands,  and  systems  of  bands  of  absorption  in  the  ultraviolet,  far  red, 
and  infrared  spectral  regions.  The  fundamental  absorption  by  the 
atmosphere  is  caused  by  the  presence  of  water  vapor,  carbon  dioxide, 
ozone,  etc.,  all  of  which  have  powerful  absorption  bands. 

Water  vapor  is  the  most  effective  absorbing  component;  Its  content 
in  the  atmosphere  can  be  sharply  varied,  depending  upon  a  whole  series 
of  factors  3uch  as  temperature,  pressure,  etc.,  [1].  The  major  quantity 
of  water  vapor  is  concentrated  at  altitudes  up  to  5  km.  Therefore  when 
determining  the  transparency  of  an  atmospheric  aerosol  it  is  necessary 
to  know  all  the  spectral  intervals  in  which  absorption  of  radiation 
by  water  vapor  is  essential.  Table  1  shows  the  basic  absorption  bands 
of  watqr  vapor  in  the  near  infrared  region  of  the  spectrum.  The 
strongest  absorption  bands  are  centered  at  the  wavelengths  1.38,  1.87, 
2.7,  and  6 . 3  ym. 

<• 

Another  absorbing  component  of  the  atmosphere  is  carbon  dioxide, 
which  is  usually  found  in  the  atmosphere  in  smaller  quantities  than 
water  vapor.  Carbon  dioxide  possesses  several  absorption  bands  in  the 
infrared  region  of  the  spectrum,  located  in  the  wavelength  region 
X  >  2,ym. 

Intensive  absorption  bands  in  the  infrared  region  are  caused  by 
vibration  transitions,  where  only  vibrations  with  frequencies  of 
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-  673  cm”1  and  v2  «  2378  cm"1  are  active  in  absorption.  The 
indicated  basic  frequencies  reflect  two  intensive  absorption  bands  — 
near  4.3  and  14. 7  ym.  Besides  these  basic  bands,  C02  absorption  bands 
centered  at  wavelengths  2.7  and  10  y  and  a  number  of  other  weak  bands 
are  also  detected.  Table  1  gives  all  significant  absorption  bands  for 
CO- 

C  t 


Table  1  , 


Absorption  bend*  | 

Transparency 
window,  pul 

0,95-1,06 

watar  vapor 

1  CO, 

designation 
of  band 

spsetral 

region, 

ym 

band 

center, 

ym 

band 

center 

width  of 
band,  em“* 

a 

0,70-0,74 

0,718 

1.4 

6650-7250 

1,2— 1,3 

'  0.8  ym 

0,79—0,84 

0,810 

1,6 

6000-6550 

M-l.8 

0,926-0.978 

0,935 

2,0 

4750-5200 

2, 1-2, 4  ' 

0  ‘ 

1.095^-1,165 

1,130 

2.7 

3480-3800 

3, 3-4.2" 

+ 

Ml  9- 1,498 

1,395 

4.3 

2160-2500 

4,5— 5,0 

u 

1,762-1,977 

1,870 

4.8 

1980-2100 

8-13 

l 

2,520—2,845 

.2,680 

5,2 

1870-1980 

3.2  ym 

2,990—3,570 

ia 

14,7- 

550-800 

•  * 

6.3  ym 

4,800—8,690 

6,3 

During  measurements  of  the  spectral  transparency  of  an  atmospheric 
aerosol  it  is  also  necessary  to  take  into  account  absorption  of  short¬ 
wave  radiation  by  ozone,  whose  concentration  shows  a  rather  unique 
variation  with  altitute:  it  increases  sharply  at  altitudes  of  22-27  km 
and  is  reduced  virtually  to  zero  at  altitudes  above  40  km.  Observations 
carried  out  in  natural  and  laboratory  conditions  have  shown  that  ozone 
has  several  absorption  bands  in  the  solar  radiation  region  of  the 
spectrum,  with  the  most  Intensive  absorption  band  in  the  ultraviolet 

O 

section  of  the  spectrum,  i.e.,  2200-3200  A.  In  the  Infrared  region  of 
the  spectrum  under  atmospheric  conditions  only  a  narrow  but  Intensive 
band  at  9*6  ym  is  observed  clearly.  The  remaining  ozone  absorption 
bands  are  overlapped  by  more  Intensive  absorption  bands  of  water  vapor 
and  carbon  dioxide.  All  of  the  ozone  absorption  bands  are  manifested 
clearly  only  under  the  condition  that  the  source  of  radiation  be 
located  outside  the  terrestrial  atmosphere  (Sun,  Moon,  stars).  If  the 
source  of  radiation  is  located  in. layers  of  the  atmosphere  below 
15  km,  absorption  by  ozone  will  be  very  Insignificant.  Consequently, 
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during  determination  of  the  transparency  of  an  atmospheric  aerosol  with 
an  artificial  radiation  source  the  influence  of  ozone  absorption  can 
be  ignored. 

Thus,  the  considered  absorption  bands  create  the  so-called  picture 
of  absorption  of  radiation  by  the  atmosphere,  a  picture  which  changes 
depending  upon  the  concentration  and  state  of  the  individual  absorbing 
components.  In  designing  optical  and  electronic  instruments  it  is 
extremely  important  to  note  the  presence  of  "windows"  of  atmospheric 
transmission  (see  Table  1). 

The  changes  which  occur  in  the  concentrations  of  the  absorbing 
components  of  the  atmosphere  make  it  necessary  to  modify  the  methods 
used  to  calculate  the  absorption  coefficient  during  distribution  of 
radiation  in  nonhorizontal  directions.  Such  methods  have  been 
developed  for  several  cases  which  are  important  in  the  practical 
respect.  For  example,  the  method  of  effective  absorbing  mass  Is  applied 
for  the  case  of  a  beam  of  pressure  which  is  constant  in  direction  [2], 

Scattering  of  radiation  arises  in  those-  cases  when  the  medium  is 
nonuniform,  i.e.,  scattering  occurs  in  the  presence  of  spatial 
inhomogeneity  of  the  dielectric  constant.  In  atmospheric  conditions 
the  basic  factors  which  determine  scattering  are  fluctuations  in  the 
density  of  the  air  and  the  presence  of  aerosol  particles  (drops  of 
water,  dust  particles,  etc.).  According  to  the  work  by  Rayleigh  [3], 
for  thecoefficient  of  molecular  scattering  we  will  have 

ap  2 V|(mi  —  1  )S  (2) 

where  mi  is  the  refractive  index  for  the  i-th  gas  under  normal 
conditions,  N  is  the  number  of  molecules  per  unit  volume,  v^  Is  the 
partial  volume  of  the  i-th  gas  in  air,  and  p.^  is  the  coefficient  of 
depolarization  of  the  given  gas.  The  values  of  the  coefficients  of 
molecular  scattering  for  Individual  wavelengths  were  calculated  by 
R.  Penndorf  [4]. 


Analysis  of  the  data  given  in  [*1]  shows  that  molecular  scattering 
can,  for  practical  purposes,  be  ignored  in  the  infrared  region  of  the 
spectrum.  This  factor  must  be  taken  into  account  in  the  visible  region 
of  the  optical  spectrum. 

However,  molecular  scattering  is  not  the  principal  factor  in  the 
scattering  of  radiation  in  the  atmosphere.  Here  there  is  always  a 
process  connected  with  scattering  on  the  aerosol  component  of  the 
atmosphere  —  the  most  variable  portion  of  the  atmosphere.  Thus, 
according  to  the  works  by  S.  F.  Rodionov  [53,  T.  P.  Toropova  [6],  and 
others,  scattering  on  aerosols  is  comparable  with  the  magnitude  of 
molecular  scattering  of  radiation  even  in  the  absence  of  fog  and 
noticeable  smoke.  In  a  turbid  atmosphere,  i.e.,  in  the  presence  of 
smoke,  scattering  on  aerosols  is  the  principal  factor  operating  to 
attenuate  radiation,  so  that  molecular  scattering  is  negligibly  small 
in  comparison  with  it.  Foltzlk  [73  investigated  the  dependence  of  the 
polydisperse  coefficient  of  scattering  on  the  size  and  number  of 
particles  of  a  dispersed  system.  He  found  that  this  relationship 
carries  a  spectral  character.  Thus,  all  investigations  which  have  been 
carried  out  attest  to  the  fact  that  the  transparency  method  can  be 
applied  to  determine  the  dimensions  of  particles  of  a  natural  aerosol. 
From  work  [83  it  is  clear  that  the  transparency  method  is  suitable  for 
application  to  disperse  systems  of  particles  with  r  <  5  pm  (in  real 
conditions  this  would  be  atmospheric  dust  and  smoke).  The  limited 
nature  of  application  of  the  transparency  method  to  disperse  systems 
of  smoke  and  dust  is  explained  by  the  fact  that  determination  of 
dimensions  of  particles  with  r  >  5  pm  requires  knowledge  of  the 
polydisperse  coefficient  of  scattering  in  a  comparatively  far  infrared 

-A 

regjJon.  Working  with  such  a  spectral  interval  in  the  real  atmosphere 
is  connected  with  major  difficulties,  both  in  principle  and  from  the 
technical  point  of  view. 

Measurements  of  the  Transparency  of  an  Atmospheric  Aerosol 

Any  measurement  of  the  transparency  of  a  certain  layer  of  the 
atmosphere  in  the  horizontal  direction  normally  reduces  to  comparison 
of  the  intensity  of  radiation  which  has  passed  through  the  layer  with 
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the  Intensity  of  the  incident  radiation  or,  in  general,  with  the 
radiation  of  a  certain  source,  yielding  a  comparison  spectrum.  The 
complexity  of  experiments  of  this  type  for  the  visible,  near  ultraviolet, 
and  infrared  regions  is  determined  by  the  necessity  of  utilizing  long 
paths,  traversed  by  the  beam,  because  of  the  great  transparency  of  air. 
Selection  of  the  magnitude  of  the  working  route  depends  essentially  on  the 
sensitivity  of  the  receiver,  since  it  is  necessary  to  deal  with 
comparatively  small  changes  in  the  intensity  of  radiation  caused  by 
scattering  on  the  natural  aerosol.  During  study  of  the  dynamics  of 
processes  occurlng  in  aerosols  the  working  path  is  selected  as  a 
constant,  usually  equal  to  5-10  km.  It  is  desirable  that  the  path  be 
at  as  great  a  distance  as  possible  from  the  surface  of  the  Earth,  in 
order  to  eliminate  random  changes  in  the  optical  state  of  the  route. 

Determination  of  transparency  g*(v*)  can  be  carried  out  as 
follows.  A  flux  of  radiation  from  a  source  located  in  the  focus  of  the 
sending  mirror  passes  through  a  layer  of  the  atmosphere  and  upon  being 
reflected  from  the  return  mirror  arrives  at  the  receiving  mirror, 
which  directs  the  flow  of  radiation  to  the  input  slot  of  a  monochromator. 
After  dispersion  the  monochromatic  beam  is  directed  to  the  sensing 
surface  of  the  receiver.  The  amplified  signal  of  a  wide-band  amplifier 
is  detected  and  is  measured  by  a  recording  device.  Another  light  flux 
(control  beam)  which  has  not  been  subjected  to  attenuation  oh  aerosols 
is  directed  by  a  system  of  concave  spherical  mirrors  from  the  same 
source  onto  the  input  slit  of  the  monochromator.  The  ratio  of  intensity 
of-the  first  beam  to  the  second  is  proportional  to  the  transmission. 

An  advantage  of  this  circuit  Is  the  fact  that  it  makes  it  possible  to 
arrange  the  source  of  radiation  close  to  the  receiving  device  and  to 
use  the  given  source  of  radiation  for  the  control  signal. 

The  intensity  of  radiation  which  has  passed  through  an  atmosphere 
containing  an  aerosol  is  determined  by  the  formula 

/.-/o*-*',  (3) 

where  1^  is  the  intensity  of  radiation  after  passage  along  path  1  and 
Iq  is  the  intensity  of  the  same  source  In  direction  l. 
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The  intensity  of  the  control  signal  equals 

(4) 

The  distance  «  l  is  selected  in  such  a  manner  that  scattering 
on  the  aerosol  can  be  considered  to  be  absent.  Taking  into  account  the 
loss  of  intensity  of  radiation  during  reflection  from  the  mirror*  we 
will  rewrite  (3)  and  (4)  in  the  form 

Ix-rM-Qi,  (5) 

/j  ■■  r'7 0*-^ ,  ( 6 ) 

where  r  and  r1  characterize  the  loss  of  intensity  of  radiation  on  the 
n  n 

path*-  l  and  l ^  due  to  reflection. 


Solving  (5)  and  (6)  simultaneously,  we  obtain 


(7) 


Here  a  preliminary  study  of  the  dependence  of  rn  and  r^  on 
wavelength  is  indicated.  Thus  the  given  method  makes  it  possible  to 
determine  the  magnitude  of  attenuation  of  radiation  with  an  accuracy 
up  to  a  certain  constant  A. 

If  the  measurement  of  intensity  is  conducted  at  the  terminal  and 
intermediate  points,  it  will  be  possible  to  determine  the  scattering 
of  radiation  by  an  aerosol  which  is  located  between  these  points.  Here 
one  should  take  into  account  the  fact  that  if  a  source  is  a  point  source 
the  reduction  in  measurable  intensity  of  radiation  with  distance  from 
the  source  occurs  not  only  as  the  result  of  absorption  and  scattering, 
but  basically  is  a  consequence  of  a  drop  in  intensity  by  the  law  of 
inverse  squares.  Calculation  of  this  type  of  attenuation  requires  the. 
application  of  neutral  filters  which  will  compensate  attenuation  due  to 
distance  [9]. 
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Measurements  conducted  in  this  way  make  it  possible  to  determine 
attenuation  of  radiation  Q  on  a  section  of  the  path  from  the  intermediate 
to  the  end  points,  i.e., 

Q-Ap  +  tfV).  (8) 

By  calculating  the  values  of  kp  by  the  Rayleigh  formula  (2),  we 
will  be  able  to  determine  the  value  of  the  polydisperse  coefficient  of 
scattering. 


If  the  source  of  radiation  is  the  Sun  the.  coefficient  of  attenuation 
is  determined  according  to  the  formula* 

q  —  In  /*,  —  In  /«,  (9) 

mt, — mt% 


where  m_  -  m  represents  the  difference  of  atmospheric  masses  between 
z2  Z1 

the  ob  ervation  points;  I_  and  I  are  the  monochromatic  fluxes  of 

Z1  z2 

direct  solar  radiation  at  the  two  altitudes  z1  and  z2  for  a  given 

zenith  distance.  The  difference  in  atmospheric  masses,  m  -  m  ,’is 

Z1  z2 

determined  according  to  the  Bemporad  table  with  a  correction  for 
pressure . 


Thus,  by  applying  one  of  the  methods  outlined  above  we  are  able 
to  determine  the  spectral  path  of  the  polydisperse  coefficient  of 
scattering,  knowledge  of  which  Is  necessary  for  calculation  of  the 
spectrum  of  dimension  of  particles  of  a  dispersed  system  by  the 
transparency  method. 


The  examined  methods  of  determining  the  polydisperse  coefficient 
of  scattering  g* (v*)  assumed  that  the  spectral  measurements  are 
conducted  outside  the  bands  of  absorption  by  the  gaseous  components 
of  the  atmosphere.  In  the  opposite  case  it  is  necessary  to  introduce 
a  correction  for  absorption,  and  since  water  vapor  is  the  most  variable 
and  the  strongest  absorbing  component,  in  the  course  of  the  experiment  it 
is  necessary  that  data  on  absolute  humidity  be  available. 
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During  measurement  of  the:quantity  g*(v*)  it  should  be  remembered 
that  formula  (3)  is  valid  for  attenuation  of  intensity  in  a  direction 
directly  ahead  of  1^.  However,  radiation  which  has  been  scattered  from 
the  side  will  be  incident  on  the  receiving  area;  this  is  connected  with 
the  angular  dimensions  of  the  photometer.  We  shall  consider  that  only 
the  role  of  once-scattered  radiation  is  great.  We  3hall  designate  as 
0  the  cone  angle  of  radiation  of  the  source,  and  \p  will  be  used  to 
designate  the  aperture  angle  of  the  receiver.  Then  for  the  intensity 
of  single  [scattered]  radiation  we  can  write 

/,  -  JJ  -p-p  *-«<;-+'*>p(»)cos  <|>  dv ,  ( 10 ) 

v  *  1 

where  P(ft)  is  the  indicatrix  of  scattering  and  r^  and  are  the 
distances  from  the  elementary  scattering  volume  to  the  source  and  to 
the  receiver,  respectively. 

For  the  total  intensity  recorded  by  the  receiver  we  have 

(11) 

C  -  -i-  f  ?  ^  +  ®)  ccs  <!> d  <|>  d  ft'. 

From  (11)  it  follows  that  for  specific  experimental  conditions 
Q  and  C  can  be  obtained  by  relative  measurements  of  I  for  three 
receiver  -  source  distances. 

Device  for  Automatic  Computation 
of  the  Distribution  Function 

In  practice  cases  are  encountered  when  it  is  necessary  to  know 
the  distribution  of  particles  with  respect  to  dimensions  immediately 
after  measurement  of  the  spectral  transparency  of  a  disperse  system. • 
Calculation  of  the  curve  of  particle  distribution  with  respect  to 
dimensions,  f(r),  can  be  significantly  accelerated  if  we  introduce 
automation  of  the  computation  process  according  to  the  equations 
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(12) 


x  ^g^Tfj^axj)  +  «e(a  *)+■£*■  «o,(a  t)  J. 


f(r) 


m{a) 
2*  fl’rj 


(13) 


Pig.  1.  Block  diagram  of  calculation 
device . 


This  is  accomplished  by  means  of  digital  or  other  analogous 
computer  technology.  Together  with  LITMO  [Leningrad  Institute  of 
Precision  Mechanics  and  Optics]  we  conducted  such  a  study,  as  a  result 
of  which  we  developed  a  computer  for  calculating  the  spectrum  of 
particles  by  the  transparency  method.  The  block  diagram  of  such  a 
device  for  calculating  the  function  f(r)  is  shown  on  Fig.  1.  After 
conversion  into  g(x/2)  in  block  C,  the  experimental  function  1(A)  is 
passed  to  the  converter  block  CB  either  directly  or  through  the  analog 
memory  AM.  The  kernel  u(ax)  is  passed  to  the  converter  block  from 
functional  unit  FU  simultaneously  with  g(x/2).  The  product  g(x/2)u>(ax) 
obtained  in  the  CB  passes  to  the  integrator,  which  develops  the  unknown 
function  f(r)  in  the  form  of  a  set  of  values  fixed  by  integrator  In. 
Control  unit  CU  provides  the  following:  synchronization  of  input 
g(x/2)  and  reproduction  of  w(aXj),  setting  x^  ,  and  control  of  the 
integration  cycle  duration  by  changing  the  scale  of  the  rate  of  input 
of  experimental  Information. 


Thus,  obtaining  the  spectrum  of  particle  distribution  with  respect 
to  dimensions,  f(r),  is  connected  with  multiple  use  of  the  experimental 
information.  A  special  analog  computer  realizes  relationships  (12) 
and  (13)  in  accordance  with  the  block  diagram  shown  on  Fig.  1.  By 
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Fig.  2.  Schematic  diagram  of  an  analog  computer. 

means  of  photoreceiver  PR  (Fig.  2)  intensity  I(X)  of  the  flow  of 
radiation  is  converted  to  a  voltage  which  varies  according  to  the 
argument  X  whi<-h  is  introduced  from  pickup  D. 

The  obtained  function  I ( X )  is  fed  to  logarithmic  amplifier  LA. 

As  a  second  factor,  the  quantity  IQ(X)  is  fed  to  the  LA  from  memory 
unit  MU.  The  obtained  value  g(X)  is  fed  to  converter  C,  where 
transition  to  g(x/2)  is  accomplished.  The  obtained  value  of  g(x/2)  is 
fed  to  multiplier  device  ,  to  the  second  input  of  which  the  value 
R(x,  r)  of  the  exact  solution  of  equation  (12)  is  fed  from  functional 
F^.  The  product  g(x/2)R(x,  r)  obtained  at  the  output  of  multiplier 
unit  is  integrated  with  respect  to  y,  in  device  E.  After  division 

of  a  generated  value  by  r^  the  ordinate  f(r)  of  the  unknown 

curve  of  particle  distribution  is  obtained  in  unit  M.  The  magnitude 
of  this  ordinate  is  recorded  for  each  discrete  value  of  r  by  indicator 
Ind.  As  the  result  of  multiple  repetition  of  the  described  process 
for  each  value  of  r  in  the  range  of  particle  dimensions  of  interest, 
the  point  curve  f(r)  characterizing  the  distribution  of  particles  of 
the  polydisperse  medium  under  investigation  is  obtained  on  the 
indicator.  Control  unit  CU  ensures  matching  of  the  work  of  all 
circuits  in  the  computer  and,  depending  on  the  specific  structure  of 
the  device,  can  fulfil  various  functions.  The  total  time  for  calculation 
of  the  spectrum  depends  on  the  computation  means  and  the  speed  of 
information  input;  it  does  not  exceed  10  minutes. 


FTD-HT-23-97-70 


11 


Thus ,  the  existence  of  a  high-speed  instrument  and  a  computer 
device  permits  operational  calculation  of  the  curve  of  distribution  of 
particles  with  respect  to  size. 
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